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[i] The phase scintillation index was obtained from LI GPS data collected with the Canadian High 
Arctic Ionospheric Network (CHAIN) during years of extended solar minimum 2008-2010. Phase 
scintillation occurs predominantly on the dayside in the cusp and in the nightside auroral oval. We set 
forth a probabilistic forecast method of phase scintillation in the cusp based on the arrival time of either 
solar wind corotating interaction regions (CIRs) or interplanetary coronal mass ejections (ICMEs). CIRs 
on the leading edge of high-speed streams (HSS) from coronal holes are known to cause recurrent 
geomagnetic and ionospheric disturbances that can be forecast one or several solar rotations in advance. 
Superposed epoch analysis of phase scintillation occurrence showed a sharp increase in scintillation 
occurrence just after the arrival of high-speed solar wind and a peak associated with weak to moderate 
CMEs during the solar minimum. Cumulative probability distribution functions for the phase 
scintillation occurrence in the cusp are obtained from statistical data for days before and after CIR and 
ICME arrivals. The probability curves are also specified for low and high (below and above median) 
values of various solar wind plasma parameters. The initial results are used to demonstrate a forecasting 
technique on two example periods of CIRs and ICMEs. 

Citation: Prikryl, P., P. T. Jayachandran, S. C. Mushini, and I. G. Richardson (2012), Toward the probabilistic 
forecasting of high-latitude GPS phase scintillation. Space Weather, 10, S08005, doi:10.1029/2012SW000800. 


1. Introduction 

[ 2 ] Interplanetary coronal mass ejections (ICMEs) and 
corotating interaction regions (CIRs) on the leading edge 
of high-speed streams (HSSs) from coronal holes are 
major causes of ionospheric disturbances. The strongest of 
these solar wind events result in major geomagnetic 
storms [Tsuriitani and Gonzalez, 1997; Richardson et ah, 2006; 
Tsurutani et ah, 2006; Badruddin, 2009], which have a global 
impact on the ionosphere [Buonsanto, 1999; Mannucci et ah, 
2005; Pokhotelov et ah, 2010; Pedatella and Forbes, 2011]. 
Although there are many differences between CIR and 
CME driven geomagnetic storms [Borovsky and Denton, 
2006], they both result in enhanced ionospheric irregular- 
ities causing radio wave amplitude and phase scintillation. 
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which can affect operation of satellite communication and 
navigation systems [Aarons, 1982; Basu et ah, 2002]. Scin- 
tillation at radio frequencies from VHF to L band has been 
intensely studied for geomagnetic storms at low latitudes 
[Crane, 1976; Dabas et ah, 1989; Aarons, 1982; Basu et ah, 
2005; Ray and DasGupta, 2007; Li et ah, 2006, 2008; de Abreu 
et ah, 2010], and at high latitudes [Aarons, 1982, 1997; 
Aarons et ah, 2000; Mitchell et ah, 2005]. Ionospheric scin- 
tillation degrades GPS positioning accuracy and causes 
cycle slips, which can lead to complete loss of signal phase 
lock and operational outages. The impact on signal track- 
ing performance of a GPS receiver was studied theoreti- 
cally [Hegarty et ah, 2001; Conker et ah, 2003] and practically 
[Skone and de Jong, 2000; Sreeja et ah, 2011; Strangeways et ah, 
2011 ]. 

[ 3 ] Ionospheric scintillation is difficult to predict for a 
given raypath because of complexity of ionospheric irreg- 
ularity generation, coupling of the solar wind magnetic 
and electric fields to the ionosphere, and variable con- 
vection of ionospheric plasma. Past research and devel- 
opment of methods for specification and forecasting of 
scintillation have focused on the equatorial region where, 
due to the presence of an equatorial anomaly, the most 
intense scintillation is known to occur after sunset [Basil 
et ah, 2002, 2010; Anderson et ah, 2004; de La Beaujardiere 
et ah, 2004; Redmon et ah, 2010]. However, statistical 
(empirical) models of scintillation have been developed for 
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Canadian High-Arctic Ionospheric Network (CHAIN) 



Figure 1. Canadian High Arctic Ionospheric Network (CHAIN). The corrected geomagnetic 
(CGM) latitudes 70° and 80° are superposed over the geographic grid. 


both the equatorial and high-latitude ionospheres [ Secan 
et al, 1995, 1997]. These models used scintillation data 
from past satellite experiments (Wideband, HiLat and 
Polar BEAR) that covered frequencies from VHF to L-band, 
although the modeling was based on VHF data only. More 
recently, using satellite in situ measurements of plasma 
(electron) density fluctuations, IRI ionosphere model, and 
the phase screen propagation model, a scintillation clima- 
tology model for the high-latitude ionosphere was devel- 
oped [Wernik, et al, 2007]. As the number of ground-based 
scintillation monitors of global navigational satellite sys- 
tems (GNSS) signals in the auroral and polar regions 
increases, it has become possible to construct scintillation 
climatology from accumulated scintillation data [Spogli 
et al, 2009; Prikryl et al., 2011a] and to address the ques- 
tions of forecasting and mitigation of L-band scintillation at 
high-latitudes [Decker and Pedersen, 2001; Aquino et al, 2007; 
Strangeways et al, 2011]. In this paper we present observa- 
tions of GPS phase scintillation following the arrivals of 
solar wind disturbances (ICMEs and HSS/CIRs) and dis- 
cuss possibilities of forecasting the phase scintillation 
occurrence at high latitudes. 

2. Instruments and Data 

[ 4 ] The Canadian High Arctic Ionospheric Network 
(CHAIN) [Jayaclrandran et al, 2009] consists of ten 


specialized GPS receivers and six Canadian Advanced 
Digital Ionosondes (CADIs) distributed in the auroral oval, 
cusp and the polar cap (http://chain.physics.unb.ca/chain). 
A brief description of CHAIN instruments, the scintillation 
parameters obtained and the geographic and corrected 
geomagnetic coordinates of stations can be found else- 
where [e.g., Prikryl et al, 2011a]. 

[ 5 ] The CHAIN (Figure 1) is equipped with NovAtel 
OEM4 dual frequency GPS receivers GSV 4004B [Van 
Dierendonck and Arbesser-Rastburg, 2004] with special firm- 
ware specifically configured to measure and record power 
and phase of the GPS LI signal at high sampling rate 
(50 Hz). The receiver computes ionospheric total electron 
content (TEC) using both GPS LI and L2 signals, the 
amplitude scintillation index S 4 and the phase scintillation 
index <r,|, that is used in this study. The phase scintillation 
index < 7 $ is the standard deviation of the detrended phase 
using a filter in the receiver with 0.1 Hz lower cutoff fre- 
quency. Intense scintillation is usually accompanied by 
cycle slips. A cycle slip [Horvath and Crozier, 2007] is 
defined here as a jump in differential phase TEC of more 
than or equal to 1.5 TECU in one second (1 TEC Unit 
corresponds to 10 16 electrons/m 2 ). 

[6] CHAIN is supported by radars, riometers, optical 
instruments and magnetometers of the Canadian Geospace 
Monitoring (CGSM) program [Liu, 2005], including the 
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YEAR: 2008 MONTH: 12 



DAYS 

Figure 2. Hourly averages of the solar wind velocity V (black line), flow angle (j) V (black 
crosses), the proton temperature T p (yellow), the interplanetary magnetic field magnitude |B| 
(red), solar wind density n p (green), and the mean standard deviation a Bz (purple dotted line) 
of the IMF B z are shown. Asterisks and triangles at the bottom of the figure indicate CIR and 
ICME arrivals, respectively. Phase scintillation occurrence (blue) observed from Cambridge 
Bay and Taloyoak is superposed. 


Northern Solar Terrestrial Array (NORSTAR) [Donovan 
etal., 2003]. 

3. Solar Wind Data 

[ 7 ] Solar wind origins of ionospheric disturbances 
causing the observed scintillation are identified as either 
CIR/HSSs or ICMEs. For ICMEs, we use a catalog of near- 
Earth ICMEs [Richardson and Cane, 2010] that is updated at 
http://www.srl.caltech.edu/ACE/ASC/DATA/level3. Solar 
wind plasma parameters are obtained from the Goddard 
Space Flight Center Space Physics Data Facility OMNIWeb 
(http://omniweb.gsfc.nasa.gov/) and are used to identify 
arrival times of CIR/HSSs [Prikryl et al., 2009b]. In brief, 
as a consequence of solar rotation, a high speed stream 
generates a CIR, a region between the fast and slow solar 
wind where the solar wind plasma is compressed [Pizzo, 
1978]. A stream interface is characterized by an abrupt 
and simultaneous drops in particle density and rises in 
proton temperature within [Belcher and Davis, 1971; 
Burlaga, 1974; Gosling et al, 1978], which is followed by a 
peak in magnetic field magnitude [Jian et al., 2006]. The 
standard deviation a Bz of the IMF B z component also peaks 
following the stream interface and remains elevated for a 
few days. These signatures are associated with bipolar 
east-west deflections of the solar wind flow [Gosling et al, 
1978]. In the geocentric solar ecliptic (GSE) coordinates, 
the flow angle (j> v increases from zero along -X G se and is 
positive toward +Y GSE . The flow angle 0 V is often used as 
the solar wind parameter to determine arrival times of CIR 
stream interfaces [McPherron and Siscoe, 2004]. For this 
study, a mid-point between peaks in plasma (number) 


density n p and magnetic field magnitude |B| at the leading 
edge of rising solar wind velocity V is used to identify the 
stream interface, the discontinuity separating slow and fast 
solar wind [Tsurutani et al, 2006]. This then provides the 
key (zero) time for superposed epoch analysis (Section 5). 
The flow signatures, which may differ from the latter by a 
few hours, were used only when peaks in |B| or n p were 
ill-defined or when these data were missing. A total of 
64 HSS/CIRs with a maximum velocity Vmax exceeding 
500 km/s were identified in the period from 2008 to 2010. 
To characterize the solar wind by one index, we use the 
quasi-invariant (QI) index [Osherovich et al, 1999] that is 
defined as the ratio of the solar wind magnetic to ram 
pressures. 

[8] Figure 2 shows solar wind plasma parameters in 
December 2008 when three major CIRs (December 3, 22 
and 31) at the leading edges of HSSs (Emax > 500 km/s), 
and a minor CIR/HSS on Decemberll, were observed. 
The CIRs at the leading edges of HSSs show peaks in n p 
and |B| and an increase in the mean standard deviation a Bz 
of the IMF B z component. The enhanced a Bz indicates 
enhanced amplitudes of the IMF fluctuations including 
solar wind Alfven waves. Also, two ICMEs arrived on 
December 4 and 16 followed by decreased proton tem- 
perature (shown in yellow lines), which is one of the 
characteristics of ICMEs [e.g., Zurbuchen and Richardson, 
2006, and references therein]. The first of these ICMEs 
closely followed arrival of the HSS a day before so it is not 
possible to distinguish its effect on scintillation from that 
of the CIR/HSS. The phase scintillation occurrence 
obtained from the Cambridge Bay and Taloyoak GPS data 
combined is shown (blue line) in the midsection of 
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CHAIN 2008-2010: SIGMA PHI >0.1 OCCURRENCE 

UIT DISTURBED DAYS n=126 
more than 60% Kp)2 


Figure 3. Phase scintillation occurrence maps for (a) low and (b) moderate Kp index. 


Figure 2. The second ICME was followed by an increase in 
phase scintillation occurrence, as were all four CIRs shown 
in Figure 2. 

4. GPS Phase Scintillation at High Latitudes 

[ 9 ] At high latitudes, the maximum occurrence of GPS 
phase scintillation was found in the auroral oval near 
magnetic midnight [Aarons, 1997] particularly during 
magnetic storms [Aarons et al, 2000]. Flowever, phase 
scintillation is known to also occur on the dayside in the 
cusp and in the polar cap [Basil et al, 1998; Spogli et al., 
2009; Alfonsi et al., 2011; Prikryl et al, 2011a, 2011c]. A cli- 
matology study for solar minimum years of 2008 and 2009 
[Prikryl et al., 2011a] showed that phase scintillation, as a 
function of magnetic local time and geomagnetic latitude, 
primarily occurs on the dayside in the ionospheric foot- 
print of the cusp and on the night-side in the auroral oval. 
As expected, these scintillation regions shift in latitude in 
response to varying geomagnetic activity. The statistics 
have been extended to 2010 [Prikryl et al, 2011b] when the 
solar activity started to increase in the new solar cycle. 
Figures 3a and 3b show the maps of the percentage occur- 
rence of phase scintillation with 0 ® exceeding 0.1 radians 
as a function of magnetic local time (MLT) and CGM lati- 
tude for 970 quiet and 126 disturbed days selected using 
Kp = 2 as a threshold. Superposed on the maps are the 
boundaries of the Feldstein statistical auroral oval for quiet 
and moderately disturbed conditions [Holzworth and Meng, 
1975]. As expected, the scintillation occurrence is confined 
in a smaller and narrower area during quiet conditions. 
The scintillation area becomes larger and broader when 
the auroral oval expands during disturbed conditions. It is 
noted that the phase scintillation occurrence map in 


Figure 3b is qualitatively similar to the observed and 
modeled height-integrated electron-density irregularity 
strength from the WideBand Model (http://www.nwra. 
com/ionoscint/wbmod.html) for moderate solar activity 
level (SNN 50 to 75), which is much higher than the aver- 
age SNN during the last deep solar minimum. The simi- 
larity in the latitudinal and MLT extent of enhanced 
scintillation in the cusp and auroral oval is also found 
when comparing with the modeling results by Wemik et al. 
[2007, Figure 13]. 

[ 10 ] Dotted concentric circles in Figure 3b show latitudes 
of Cambridge Bay (CBB) and Taloyoak (TAL). The CBB 
data were available since December 2007 while TAL star- 
ted operation in July 2008. The combined CBB + TAL data 
set samples the poleward edge of auroral oval in the night 
[Prikryl et al, 2010] but is representative of the cusp iono- 
sphere during the daytime. The combined CBB + TAL 
time series of phase scintillation occurrence are used to 
examine correlation with solar wind parameters and to 
examine the possibility of forecasting phase scintillation 
following CIR/HSSs and ICMEs. 

[ 11 ] We now compare the scintillation occurrence over 
the CFLAIN array during days of active ICMEs (Figure 4a) 
and HSS/CIRs (Figure 4b) with the geomagnetically dis- 
turbed days based on Kp index as shown in Figure 3b. The 
catalog of near-Earth ICMEs [Richardson and Cane, 2010] 
includes 27 ICME events from 2008 to 2010 (50 days in 
total). It is noted that only a few of these ICMEs caused 
significant geomagnetic storms (minimum Dst of —123 nT) 
and intense scintillation events [e.g., Prikryl et al, 2011c]. 
As shown in Figure 4a, the phase scintillation during these 
50 ICME days was observed mainly in the cusp and in the 
auroral ovai. As the auroral oval expanded during dis- 
turbed conditions after the arrival of ICMEs the area of 
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CHAIN 2008-2010: SIGMA PHI >0.1 OCCURRENCE 



Figure 4. Phase scintillation occurrence maps for (a) ICMEs and (b) CIRs/HSSs. 


scintillation occurrence become broader. Also, the per- 
centage occurrence (Figure 4a) was significantly greater 
than for all quiet days (Figure 3a), particularly in the cusp, 
but it was about half of the occurrence for all geomagnet- 
ically disturbed days (Figure 3b). For 128 CIR/HSSs days 
(two days at the start of each 64 HSSs, i.e., the CIR day and 
the following day), the phase scintillation occurrence 
(Figure 4b) was higher than for ICME days (Figure 4a) 
but lower than for all geomagnetically disturbed days 
(Figure 3b), both in the cusp and in auroral oval. Clearly, 
the geomagnetically disturbed days included other scin- 
tillation active days than the major CIR/HSS days. These 
include the strongest ICMEs, other disturbed days during 
major HSSs after the CIR + 2 day, and some minor CIR/ 
HSS (Vmax < 500 km/s) like the one on December 10, 2008 
(Figure 2). 

5. The Superposed-Epoch Analysis of Solar 
Wind Plasma Parameters 

[ 12 ] As shown in Figure 2, arrivals of high-speed streams 
(HSS) originating from coronal holes are followed by sig- 
nificantly enhanced phase scintillation occurrence partic- 
ularly in the cusp, usually continuing for a few days. To 
obtain the mean response of scintillation relative to CIR 
arrival time we applied superposed epoch (SPE) analysis 
on time series of scintillation occurrence and solar wind 
parameters using key times of 64 CIRs. Figure 5a shows 
typical mean variations of the solar wind parameters over 
12 days centered at the key time (day 0). The mean solar 
wind velocity decreases to a minimum just before the key 
time and then steeply rises to a maximum two days later. 
All other mean solar wind parameters peak at, or very 


close to, the key time. The mean flow angle cf) V shows a 
bipolar east-west deflection. Figure 5b shows that the 
mean phase scintillation occurrence steeply rises at the 
CIR arrival day tapering off a few days later, which is 
similar to the mean QI and Kp indices. Also, the average 
number of cycle slips, which is nearly zero a few days 
before the CIR key day increases at the CIR interface and 
remains elevated during the HSS. A strong 1-day period- 
icity in scintillation CBB + TAL time series is due to cusp 
location of these stations where the scintillation maximizes 
and is only partly smoothed out in this relatively small 
statistical sample of 64 CIRs. 

[ 13 ] Figures 6a and 6b show the results of the superposed 
epoch (SPE) analysis for 27 ICMEs. Because the ICMEs are 
often preceded by shocks, the mean V, n p , |B| and a Bz 
increase near or just before the key time although the peak 
values are much lower than for CIRs. The QI index 
increases more gradually and to higher values than for 
CIRs. This is because many ICMEs are actually magnetic 
clouds and the ratio of the magnetic to ram pressures is 
higher inside the magnetic cloud. The Kp index peaks at a 
mean value of only 2. Also, the mean response of cycle slip 
number and phase scintillation occurrence relative to 
ICME arrivals is weaker than for CIRs for the full data set 
of ICMEs. However, if weak ICMEs are eliminated by 
selecting more geo-effective ICMEs (Dst < -10 nT), the 
increase in the mean phase scintillation occurrence is 
about 1.5 times larger (Figure 6c) than when no selection is 
made, and comparable to the increase after major CIR/ 
HSSs (Figure 5b). The mean number of cycle slips also 
increases at the ICME arrival although the reduced num- 
ber of epochs averaged is too small to suppress other 
unrelated random events. If the SPE analysis is further 
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days from CIR (E„„ > 500 km/s) 



days from CIR (V„^ > 500 km/s) 


Figure 5. Superposed epoch analysis of (a) solar wind 
plasma parameters V, n p , |B|, a Bz and d v and (b) QI, 
Kp, number of cycle slips N cys and phase scintillation 
occurrence from the Cambridge Bay and Taloyoak 
scintillation data combined. All time series are keyed 
by the times of major CIRs. 


limited to only 7 strongest ICMEs (Dst < -30 nT) (not 
shown) the peak mean scintillation occurrence and 
Kp index are ~12% and 3.5, respectively. 

6. Recurrent HSSs and Daily Occurrence 
of Phase Scintillation in the Cusp 

[ 14 ] Systematic behavior of recurrent HSSs that originate 
from coronal holes is well known and has become the basis 
of forecasting methods of space weather including geo- 
magnetic activity [McPherron and Siscoe, 2004; McPherron 
and Weygand, 2006]. The time separation of successive 
streams can often be determined by the heliospheric lon- 
gitude separation of equatorial coronal holes (heliographic 
latitude < 40° ), which is commonly 360 0 In, where n = 1, 2, 3 
or 4. There may be 1 to 4 recurrent coronal holes per solar 
rotation (~27.3 days as viewed from the Earth) and thus 
typical time separations of major recurrent plasma streams 
are ~27, 13.5, 9 and 7 days. As a result, such periodicities 
are often found in the solar corona, solar wind as well as in 
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Figure 6. Superposed epoch analysis of (a) solar wind 
plasma parameters V, n p , |B|, er Bz and d) V and (b and c) 
QI, Kp, number of cycle slips N cys and phase scintillation 
occurrence from the Cambridge Bay and Taloyoak 
scintillation data combined. The time series are keyed 
by the start times of ICMEs. 


6 of 16 


S08005 


PRIKRYL ET AL.: GPS PHASE SCINTILLATION 


S08005 


a OMNI 

SOLAR WIND VELOCITY (km/s) 


200 300 400 500 600 700 


2011 


2010 


co 

K 

< 

H 

>- 


2009 


2008 



3 6 9 12 15 18 21 24 27 

27-DAY SOLAR ROTATION 


b CHAIN: CBB + TAL 

SIGMA PHI >0.1 OCCURRENCE (%) 

0 3_ 6 9 12 


2011 


2010 


co 

« 

< 

w 

>- 


2009 


2008 



3 6 9 12 15 18 21 24 27 

27-DAY SOLAR ROTATION 


Figure 7. Daily averages of (a) solar wind velocity and (b) phase scintillation occurrence in 
Cambridge Bay and Taloyoak. 


the geophysical data including geomagnetic, ionospheric 
and sometimes even atmospheric indices and parameters 
[e.g., Mlynczak et ah, 2008; Thayer et al., 2008; Prikryl et ah, 
2009a; Sojka et ah, 2010; Emery et ah, 2011; Verkhoglyadova 
et ah, 2011; Love et ah, 2012]. However, differential solar 
rotation is well known for the photosphere, and can also 
manifest itself in the corona [see, e.g., Shnirman et ah, 2010, 
and references therein]. While the differential rotation of 
coronal holes has been confirmed in general, a different 
behavior was found between the isolated equatorial coro- 
nal holes and the extensions of polar holes. A rotation rate 
in the corona does increase toward the equator below 40° 
of heliographic latitude and can be more rigid there than 
in the photosphere but it increases toward the poles above 
the same heliographic latitude [Navarro-Peralta and 
Sanchez-Ibarra, 1994; Insley et ah, 1995]. 

[is] Recurrent HSSs and the response of the ionosphere 
can be best viewed in a stack plot of time series every 
27 days. Figures 7a and 7b show daily averages of solar 


wind velocity and phase scintillation occurrence for years 
2008 to 2010 starting on January 1, 2008. Two recurrent 
HSSs that persisted for more than a year and a half are 
seen in Figure 7a around days 6 and 16. In 2009, they were 
diminished in velocity as the equatorial coronal holes 
started to shrink, while midlatitude coronal holes formed 
[' Tsurutani et ah, 2011, and references therein]. In 2010, 
this regular structure has changed as the activity of the 
new solar cycle 24 including CMEs has begun to appear. 
However, recurrent streams from slower rotating coronal 
holes can be seen as tilted bands in the solar wind velocity 
plot. In 2010, there were some isolated coronal holes at 
high-latitude, as well as polar coronal holes that extended 
to lower heliospheric latitudes that were joined at helio- 
graphic equator. In contrast, the HSSs in 2008 originated 
mostly from large recurrent equatorial coronal holes 
resulting in high peak velocities and large amplitude 
Alfven waves [Tsurutani et ah, 2011]. 
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Figure 8. Cross- and auto-correlation of daily averages of phase scintillation occurrence and 
solar wind velocity. 


[i6] In response to solar wind HSSs the enhanced phase 
scintillation at CBB and TAL showed a very similar tem- 
poral pattern (Figure 7b). From January 2008 to May 2009, 
the scintillation occurrence was closely associated with two 
major recurrent HSSs although other minor HSSs per- 
sisted for several rotations. In 2010, scintillation was again 
associated with HSSs but this time the HSSs co-rotated 


with high-latitude coronal holes at about 28-day rotation 
period. 

[ 17 ] The overall (2008-2010) cross-correlation between 
the average daily solar wind velocity and the scintillation 
occurrence (Figure 8a) not only shows a peak of 0.46 at 
lag 0 but also peaks at multiples of solar rotation. 
Likewise, the autocorrelation functions (Figure 8b) for 
both time series show the presence of solar periodicity of 
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Figure 9a. The daily solar wind velocity over 6 solar rotations from September 28, 2008 to 
March 9, 2009 (first panel). The second, third, and fourth panels show "forecasts" by shifting 
the time series in time by 27, 54, and 81 days. The fifth panel is a superimposition of all the 
above panels. 


approximately 27 days. These peaks and two subsidiary 
peaks, i.e., a triple peak per solar rotation, are consistent 
with occasional presence, on average, of 3 active HSSs 
(coronal holes) per solar rotation (Figure 9a). The 9-day 


recurrence of geomagnetic activity correlated with solar 
wind velocity was most prominent in 2008, with the 7-day 
(more precisely 6.7-day) periodicity present in the pre- 
ceding years, particularly in 2006 [Love et ah, 2012]. 
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Figure 9b. The daily occurrence of phase scintillation over 6 solar rotations from September 
28, 2008 to March 9, 2009 (first panel). The second, third, and fourth panels show "forecasts" 
by shifting the time series in time by 27, 54, and 81 days. The fifth panel is a superimposition 
of all the above panels. 


[is] The feasibility of forecasting the solar wind condi- 
tions and scintillation occurrence for one to multiple solar 
rotations ahead during a period of stable recurrent HSSs 
streams is illustrated in Figures 9a and 9b. During this 


period two major and a minor HSSs were observed. By 
shifting the time series successively by 27 days the past 
data can be used as "forecasts" for multiple solar rotations 
ahead. Thus "forecasted" major and even some minor 
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PROBABILITY OF PHASE SCINTILLATION 



PROBABILITY OF PHASE SCINTILLATION 



Figure 10. The cumulative probability distribution 
functions for the phase scintillation occurrence com- 
bined for TAL and CBB during (a) slow and fast solar 
wind up to two days before and after CIRs (n = 64), 
and (b) two days before and after ICMEs (Dst < 10 nT; 
n = 18). The overall probabilities for each case (solid 
lines) are bracketed by a pair of curves (broken lines) 
for data corresponding to low and high (below and 
above median) solar wind parameter values. 


peaks in solar wind velocity and phase scintillation 
occurrence show reasonably good agreement with the 
observed (unshifted) time series up to 3 solar rotations 
ahead. 

7. Probabilistic Forecasting of Phase Scintillation 
in the Cusp Based on the SPE Analysis of Solar 
Wind CIRs and ICMEs 

[ 19 ] McPherron and Siscoe [2004] and McPherron and 
Weygand [2006] demonstrated a new promising technique 
to forecast geomagnetic disturbance indices relative to the 
passage of the stream interface between fast and slow solar 
wind. In particular, they showed the geoeffectiveness of 
what they called "corotating compression ridges," their 
reference to CIRs. McPherron and Siscoe [2004] construed 
their method of probabilistic forecasting of geomagnetic 
activity due to specific solar wind conditions as an analogy 
to traditional meteorology approach in weather forecast- 
ing based on atmospheric air mass climatologies. In this 
analogy, a stream interface between fast and slow solar 
wind (a CIR or an ICME magnetosheath), which leads to 
geomagnetic storms, corresponds to a boundary (weather 
front) between atmospheric air masses where the weather 
storms develop. By applying the "air mass" concept to 
space weather forecasting McPherron and Siscoe [2004] 
provided proof of concept of the method and discussed 
steps to improve it. 

[ 20 ] This method is readily applicable to other manifes- 
tations of space weather in the geospace environment that 
are sensitive to arrivals of HSSs, including ionospheric 
scintillation. Figure 10a shows the cumulative probability 
distribution function (PDF) for the phase scintillation 
occurrence at TAL and CBB contrasting the periods of 
calm [Borovsky and Steinberg, 2006] during two days of slow 
solar wind just before the major CIR/HSSs arrivals, and 
two "stormy" days of fast solar wind disturbances that 
followed. The probability that phase scintillation will 
exceed a given value as plotted on the abscissa is shown. 
The overall probability curves for both the slow and fast 
streams (shown by solid lines) are bracketed by probabil- 
ity curves (broken lines) for scintillation data subdivided 
by corresponding solar wind low and high (below and 
above median) values of V, |B| and a Bz . Likewise, Figure 10b 
shows the probabilistic curves for the SPE analysis of 
18 ICMEs (Dst < —10 nT) for two days before and after 
ICME arrivals. 

[ 21 ] As can be seen in Figure 10a, probability of scintil- 
lation occurrence is low before but significantly increased 
immediately after the CIR arrival. The probability of scin- 
tillation occurrence also depends on solar wind plasma 
parameter values. The PDFs results are very similar for 
days after ICMEs (Figure 10b), although the differences 
between the probabilities before and after ICME arrival are 
smaller. Unlike CIRs, ICMEs are not necessarily preceded 
by a particularly quiescent period. For example, there may 
be ongoing activity associated with a preceding ICME, and 
enhanced activity may arise from the interaction of an 
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Figure 11. The cumulative probability distribution 
functions for the phase scintillation occurrence com- 
bined for TAL and CBB stepped by 2-day intervals 
before and after (a) CIRs (n = 64) and (b) ICMEs 
(Dst < 10 nT; n = 18). 


ICME with the background solar wind. Figures 11a and 
lib show PDFs stepped by 2-days intervals before and 
after CIRs and ICMEs. The probabilities are clearly highest 
for two days after CIR and ICME arrivals. Daily PDFs can 
be generated but for clarity only 2-day PDFs are shown. 

[ 22 ] Because the effects of HSS/CIRs can last for several 
days it is expected that the scintillation would also persist 
longer than after the arrival of ICMEs. This and the contrast 
between the scintillation level during "calm" and "stormy" 
days before and after CIR arrivals is demonstrated in 
Figures 5b and 11a. The ICME effects can be very intense 
but are usually short-lived, with the intensity of the ICME 
impact and its duration depending on the IMF direction 
and the ICME speed. Because of a very small number of 
epochs superposed. Figures 6b and 6c are too noisy to show 
that. As a result, the PDFs shown in Figure lib do not 
reliably represent the probabilities for "quiet" days before 
ICME, the mean peak probabilities on the day of ICME 
arrival, or the actual duration of the ICME effect. Larger 
data sample of moderate to strong ICMEs is needed to 
provide that and to suppress the contamination by HSS/ 
CIRs and other disturbances. 

8. Discussion 

[ 23 ] In this paper, we applied the forecasting method by 
McPherron and Siscoe [2004] to forecasting phase scintilla- 
tion occurrence in the high-latitude ionosphere typical of a 
magnetospheric cusp footprint during a solar minimum. 
These results demonstrate the feasibility of probabilistic 
forecasting of scintillation; however, more extensive data- 
bases with improved spatial coverage in latitude/longitude 
to include polar cap and auroral oval, as well as coverage 
of a complete solar cycle, are needed. In particular, the 
majority of ICMEs during the last solar minimum caused 
only weak geomagnetic disturbances and none resulted in 
a major geomagnetic storm. Thus this sample of ICMEs is 
clearly not representative of periods of higher solar activity. 
In comparison, strong and recurrent high-speed streams 
from coronal holes remained active during most of the 
period studied in this paper. The obtained mean proba- 
bility curves for days before and after CIR arrival are 
expected to be similar throughout the solar cycle, although 
they may reflect the differences between slow and fast solar 
wind that are likely to be enhanced during the declining 
phase when the recurrent streams are dominant. 

[ 24 ] When a larger statistical sample becomes available 
the concept of probabilistic forecasting of scintillation will 
be further refined by compiling probability curves sepa- 
rately for different seasons, levels of solar activity, various 
threshold values of 0 $, as well as different magnetic lati- 
tudes and local times. It will be important to distinguish 
between the periods of southward and northward IMF. 
The IMF B z component (in the geocentric solar magneto- 
spheric coordinate system, GSM) controls the coupling to 
the magnetosphere and the energy deposition into the 
ionosphere. Considering the Russell-McPherron effect 
[Russell and McPherron, 1973], CIRs can be classified by 
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Figure 12. Examples of probabilistic forecasting of 
phase scintillation occurrence following the anticipated 
arrivals of (a) HSSs/CIRs and (b) ICMEs. The values of 
the cumulative probability distribution functions 
for ±6 days relative to the solar wind disturbance 
arrival day are used to obtain daily probabilities of 
phase scintillation occurrence to exceed 10%. The 
observed hourly occurrence of phase scintillation is 
shown in the bottom plots. 


IMF polarity as shown by Miyoshi and Kataoka [2008] and 
adopted by Prikryl et al. [2009a, 2009b]. For months around 
the spring and fall equinoxes, the mean IMF B z component 
systematically shifts southward near the CIR interface for 
the group of events classified as type SBZ (southward B z ), 
while the CIRs are associated with a northward shift in B z 
for the group classified as NBZ (northward B z ). If adopted 
for probabilistic forecasting of scintillation, this classifica- 
tion of CIRs should result in higher and lower probabilities 
for SBZ and NBZ groups, respectively. 


[ 25 ] Although in this paper we focus on high latitudes, it 
is appropriate to comment on the HSS effects at the low 
latitudes. It has been shown that TEC at low and midlati- 
tudes is affected by HSSs [e.g., Verkhoglyadova et al, 2011]. 
During disturbed periods the Joule heating at high lati- 
tudes controls neutral wind circulation that impacts the 
ionosphere at all latitudes and affects storms development. 
Also, there is direct penetration of the electric field to low 
latitudes as well as equatorward propagating traveling 
ionospheric disturbances (TIDs), both of which can result 
in enhanced ionospheric irregularity production at low 
latitudes. Although the occurrence of post-sunset scintil- 
lation at low latitudes is known to be more unpredictable 
than scintillation at high latitudes, the proposed method 
may have some value in forecasting of scintillation at low 
latitudes. This could be investigated by applying the 
SPE analysis on low-latitude scintillation time series. 

[ 26 ] Probabilistic forecasting of scintillation relative to 
solar wind disturbances ultimately requires accurate pre- 
diction of CIR or CME arrivals. Space weather forecasting 
has benefited from recent advances in MHD modeling of 
CMEs and CIRs [Odstrcil et al, 2005; Wu et al, 2011] as well 
as multispacecraft tracking of CIRs propagation [Williams 
et al, 2011; Davis et al, 2012]. Combining the solar wind 
modeling and monitoring efforts with improved space 
weather forecasting will lead to more realistic forecast of 
ionospheric response including scintillation occurrence 
and other stochastic phenomena modulated by solar wind. 

[ 27 ] To illustrate a possible application of PDFs shown in 
Figures 10a and 11a in forecasting scintillation occurrence 
we choose three CIRs identified in hourly OMNIWeb data 
in 2008 on September 30 (~17 UT), October 11 (~8 UT) and 
October 19 (22 UT) from the interval shown in Figure 9. As 
already discussed in Section 6, such stable recurrent HSSs 
streams can be forecast at least for one solar rotation, 
particularly in the absence of ICMEs. Applying the 
method of probabilistic forecasting one can propagate 
daily probabilities of scintillation occurrence to exceed, 
for example, 10% for ±6 days relative to October 27, 
November 7 and 15 (0 UT). This would provide a scintilla- 
tion occurrence forecast for a 27 day period composed of 
forecasts for individual CIRs (Figure 12a). For days with 
multiple PDF values due to overlapping 12-day intervals 
the maxima may be taken as a forecast for a given day. 
Individual forecasts could be updated and fine-tuned at a 
later date as more data become available, as coronal holes 
reappear on the east solar limb, solar wind models are 
updated and eventually the CIRs are observed upstream of 
the solar wind (e.g., ACE). The observed hourly scintilla- 
tion occurrence is shown in Figure 12a (bottom). The 
corresponding daily scintillation occurrence values that are 
shown in Figure 9b (first panel; days 298-325) are consis- 
tent with daily probabilities expressed in percentages. 
Clearly, the simple approach using 27-day recurrence 
period does not produce quite satisfactory result for the 
first CIR, in which case the scintillation is predicted one 
day too early. This could be improved if more precise time 
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of CIR is considered but variability in the rotation period of 
coronal holes is expected, as discussed in Section 6. 

[ 28 ] In principle, similar forecasts of scintillation proba- 
bilities for ICMEs could be attempted and added to CIR 
forecasts, of course, with shorter lead times. As an exam- 
ple, we consider recent events of two CMEs triggered by 
X-class solar flares in 2012 on March 5 at 03:30 UT and on 
March 7 at 00:02 UT. The solar wind including the ICMEs 
was modeled at the Goddard Space Weather Lab http:// 
iswa.gsfc.nasa.gov/iswa/iSWA.html. The model predicted 
the impact of the ICMEs at the Earth for early UT hours on 
March 7 and 8, respectively (http://earthsky.org/space/ 
another-major-solar-flare-during-night-of-march-6-7-2012). 
Using the CELLAS/MTOF proton monitor on the SOHO 
spacecraft (http://umtof.umd.edu/pm/), the approximate 
travel times of the first and second ICME were ~49 and 
35 h, respectively. Based on solar wind data projected to the 
bow shock (1-min OMNIWeb data), the magnetospheric 
impacts occurred on March 7, 04:30 UT and March 8, 
11:00 UT. However, in this exercise of nowcasting we adopt 
the solar wind model predictions of the ICME impact times 
on March 7 and 8, consult the statistical probabilistic curves 
(Figure lib) to obtain daily probabilities of interest and use 
them to "predict" daily occurrence of phase scintillation in 
the cusp. Figure 12b shows two identical sets of daily 
probabilities of scintillation occurrence to exceed 10% 
for ±6 days centered on March 7 and 8 superposed in blue 
and red colors, respectively. The observed hourly scintil- 
lation occurrence at Cambridge Bay is shown in Figure 12b 
(bottom) (no data were available from Taloyoak). It should 
be noted that the probabilities in this example are based on 
a very small sample of ICMEs (Figure lib), in which very 
few ICMEs matched the magnitude of the March 2012 
ICMEs. More realistic forecasts will require extensive 
database covering at least one solar cycle, which may pro- 
vide a possibility to discriminate among the intensity levels 
of CMEs. 

9. Summary and Conclusions 

[ 29 ] Based on the climatology of GPS phase scintillation 
at high latitudes that showed the scintillation predomi- 
nantly occurring in the ionospheric footprint of the cusp, 
a probabilistic forecasting method focusing on the cusp is 
proposed. The GPS data from two cusp stations in Cam- 
bridge Bay and Taloyoak in Canada are used to construct 
time series of scintillation occurrence that is correlated 
with solar wind velocity. Using a superposed epoch analy- 
sis it is shown that there is a significant increase in phase 
scintillation and cycle slip occurrence just after the arrival 
of high-speed solar wind CIRs and ICMEs. The mean 
scintillation activity in the cusp tapers off a few days later. 
Cumulative probability distribution functions for the 
phase scintillation occurrence in the cusp during solar 
minimum are obtained for days before and after CIR and 
ICME arrivals. The probability curves can also be specified 
for low and high (below and above median) values of var- 
ious solar wind plasma parameters. These initial results are 


used to demonstrate a technique to forecast phase scintil- 
lation occurrence in the cusp relative to stream interface 
between fast and slow solar wind and/or shocks associated 
with CIRs and ICMEs. A larger statistical sample is needed 
to refine the probabilistic forecasting for different seasons, 
levels of solar activity, solar wind velocity magnitude and 
the IMF orientation, different scintillation levels, and dif- 
ferent magnetic latitudes and local times. The method 
should be applicable to auroral zone and this will be fur- 
ther investigated as the solar activity increases to solar 
maximum and the coverage of auroral latitudes is 
improved in the coming years. More work on the validation 
of forecasting based on cumulative probabilistic curves of 
phase scintillation occurrence is required before the 
method can be adapted and included as a sub-product in 
space weather prediction procedures. 
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